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I. INTRODUCTION 

Military diesel engines are widely utilized in vehicles and for mobile and 
stationary power generating plants. Uncertainties in fuel resources and 
variability in fuel costs require that future diesel engines operate with 
improved efficiencies while maintainina low smoke emissions with fuel that may 
be low in quality. Thus, a better understanding of fuel injection and 
combustion processes in diesel combustion chambers is needed and strategies to 
improve performance are desired. 

Combustion in diesel engines is significantly affected by fuel spray 
injection and flow characteristics (Heinein 1 ). Many methods for improved mixing 
and dispersion of the spray to reduce local fuel rich regions have been 
investigated for improved combustion and reduced soot formation. In general, 
the breakup of fuel jets in diesel engines is accomplished by high pressure 
injection In the range 20 MPa, 200 atm (for indirect injection) up to about 140 
MPa, 1400 atm (for direct injection). The mechanism of breakup of a fuel jet at 
high velocity injection has been investigated and reviewed by Reitz and Bracco 2 
and Lin and Kang 3 . Swirl is imparted to the flow field in the chamber to 
further assist the fuel-air mixing process, and soot formation is strongly 
dependent upon swirl characteristics. However, fuel -rich regions apparently 
still persist in current engines, and methods for improved mixing and reduced 

soot formation continue to be of interest. 

Techniques for improving in-cylinder combustion characteristics include 
mechanical, chemical, and electrical methods. Mechanical methods include the 
adjustment of injection timing and changes in injection geometries (Pischinger 
and Cartellieri 4 ; Khan et al 5 ). Chemical concepts usually rely upon fuel 
additives or the application of catalysts on cylinder surfaces (Gaffney et al 6 ). 
Recently, studies of microexplosions in fuel -water and fuel -alcohol mixtures 
have led to the possibility that such a breakup process could be utilized to 
enhance dispersion of the injected fuel (Wang and Law^). Electrical concepts 
have been based on observed electrical phenomena in the combustion process. 
Electrical fields have been applied to the combustion chamber (Ahmad et al 6 ; 
Miller et al 9 ) however, the fundamental phenomena was not studied. Also plasma 
devices have been proposed to enhance the combustion process. 

The concept studied in this investigation involves the injection of a fuel 
jet that is electrostatically broken up and dispersion enhanced by the repulsion 
of like-charged droplets. Numerous attempts have been made to develop 
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electrostatic fuel atomizers for use in combustion systems and also to develop 
models to explain the electrostatic processes of breakup and dispersion of fuel 
jets, and those investigations were reviewed by Kelly* 0 . However, electrostatic 
dispersion of a fuel spray has not proven practical due to the very high 
resistivities of liquid hydrocarbon fuels (10 10 - 10 12 ohm-m) (Kelly 10 ’ 11 ; 
Luther 12 ) which presents difficulties in producing electrostatically charged 
sprays at flow rates above 0.1 m i/s (Luther 12 ; Newgard and Noon 13 ). The 
electrostatic spray device conceived and demonstrated by Kelly 10 ’ 11 has overcome 
many of these previous limitations. This Spray Triode can charge mineral oil 
Marcol -87 (viscosity - 25 cp) up to 1.8 C/m 3 for flow rates up to 1.25 mi/s 
(Kelly 10 ). 

The objectives of this research were to investigate experimentally the 
characteristics of charged diesel fuel jet atomization, dispersion and 
combustion, and to theoretically model dispersion and evaporation of charged 
fuel sprays. 

II. RESULTS 

The results of this investigation are reported in a series of papers Refs. 
A-E listed in the Table of Contents. Refs. B and E are included as Appendices 
in this report since they have not been published yet. These papers are briefly 
described and the important results are summarized as follows. 

Experimental 

Ref. A presents results obtained at one atmospheric pressure and room 
temperature. In this work a Spray Triode with two orifice sizes, 173 and 422 pm 
was tested at steady flow rates of 0.25-1.0 mi/s with mineral oil (viscosity » 
25 cp) under non-combusting conditions, and with no. 2 diesel fuel (viscosity - 
6 cp) under combusting and non-combusting conditions. The experimental setups 
for the noncombusting and combusting experiments are shown in Fig. 5 of Ref, A, 
along with a detailed diagram of the Spray Triode, Fig. 6. Spray charge 
densities were determined by measuring the current carried by the fluid with an 
electrometer in series between the collector electrode and ground. Charge 
density was determined by dividing the electric current by the measured liquid 
volume flow rate. Characterization of the spray was carried out by flash 
illumination and photography. Analysis of the spray pictures was carried out 
using an image processing system. 
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The results from Ref. A are summarized as follows: 

1. Continuous, stable fuel sprays were generated with charge densities of 1.5 
- 2.0 C/m 3 of fluid. 

2. Calculations using a model proposed by Bellan 14 showed that such charge 
densities may enhance spray dispersion under diesel engine conditions. 

3. Measured mean drop diameters were near 150 /im, being between Kelly's 

statistical equilibrium end state model 1 ^* 1 ^ and the Rayleigh limit. 

4. The electrical power required to generate these sprays was less than 10' 6 
times the chemical energy available from the fuel. 

5. The results showed considerable differences in spray characteristics 

between a diesel injector (alse studied) and electrostatic injection. 

6. Ignition and stable combustion of electrostatically dispersed diesel fuel 

jets were achieved. 

Ref. B (Appendix I) presents recent results obtained at various back 

pressures to 600 psig using a spray triode designed by the Parker Hannifin 
Company to fit into a small Volkswagon diesel engine. This spray triode has a 
254 **ni diameter orifice and can be operated up to a fuel pressure of 3000 psig, 
electric potentials of -20 kV and temperatures of 200°C. The experiments were 
conducted in a high pressure chamber with a relatively large window to allow 
spray development to be observed (Fig. 1, Ref. B). Provision was also made to 
measure the total current to the cathode if and the collector current if (charge 
carried by the fuel per unit time) by using a switching arrangement. The 
electric charge density of the diesel spray was measured over a range of flow 
rates from 0.5 to 2.5 mi/s at various chamber pressures of 0, 300 and 600 Dsig 
by applying different electric potentials to -18 kV. Characterization of the 
spray was carried out by flash illumination and photographs with and without 
applied electric fields at fuel injection velocities of 20 and 40 m/s, and at 
various voltages. Droplet sizing by image processing was done for the data at a 
chamber pressure of 300 psig. 

The results from Ref. B are summarized as follows: 

1. The average charge density of fuel drops linearly increased up to 1.5 C/m 3 
with increasing electric potential until breakdown occurred. 

2. After breakdown the charge density was reduced by 40 to 60% and again 
increased with a smaller slope as the electric potential was increased. 

3. At higher flow rates breakdown occurred at higher voltages. 
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4. At high back pressure lower charge dsnsity was obtained (about 1 C/m® at 
600 psig) and breakdown occurred at higher voltages. 

5. In general, with electric charge rapid breakup of the jet was achieved 
without large droplets in the core region. 

6. The spray angle increased to about 180 deg. with increase of electric 
potential which enhanceo dispersion of the fuel near the nozzle exit and 
reduced the dense core region. 

7. It was observed that the size of drops may be determined by the charge 
density rather than aerodynamic forces. 

8. The average drop size was about the same as the spray triode orifice 
diameter and was between the Kelly theory and the Rayleigh limit. 


In conclusion, these experimental investigations show that electrostatic fuel 
injection can be achieved at practical flow rates, and that significant early 
breakup and dispersion of charged fuel jets occur at high back pressures 
compared to aerodynamic breakup ana dispersion. These results are encouraging 
for engine injection, but demonstration of this electrostatic technique under 
pulsed fuel injection and autoignition conditions requires subsequent 

‘investigations to evaluate the potential for improvement in combustion and 
reduced soot formation. Since we haye experienced recurrent failure of the 

insulation in the high pressure electrostatic spray triode, a more durable unit 
is also needed. 


Theoretical Modeling 

In Ref. C and 0 a model of convective spray evaporation has been developed 
for dense and dilute clusters of drops. The model is for a cluster of 

monodisperse, single component, spherical drops that is exposed to a surrounding 
convective gas flow. The model takes into account droplet interactions within 
the cluster and the resulting limitations on evaporation as fuel vapor 
accumulates within the cluster. Results obtained from this model appear in Ref. 
C and 0 and, are not summarized here. This work was jointly supported by the 
Army Research Office, the U. S. Dept of Energy/Energy Conversion and Utilization 
Technology Program, and the Air Force Office of Scientific Research/Directorate 
of Aerospace Sciences. 

The model formulation for the evaporation of a cluster of monodisperse fuel 
drops was extended more recently to include the effects of drop charge. Two 
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different geometries were considered: (1) a spherical cluster of drops traveling 
in an ambient gas, and evolving in time, and (2) a steady ax i symmetric jet 
flow where axial distance corresponds to time (Lagrangian coordinate.:) 'i..- 
momentum equations for drops and gases were generalized to account for d ;; -'tip 
and radial cluster expansion. Combined with the momentum equations, the c •« < r 
energy equation becomes the determinant for cluster radius. When the dr-u r 
Rayleigh limit is reached (drop instability due to charge), th° calculation is 
stopped since the remaining mass of the drops is negligible for the sizes of 
drops of interest. 

Results obtained for the evaporation and dispersion of a spherical cluster 
of charged drops moving in an ambient gas are given in Ref. E. (Appendix II) 
For dense clusters of drops, it was found that electrostatic effects dominated, 
and that turbulence effects, modeled as described in Ref. D, were not the 
controlling parameter as was the case for uncharged drops. The charge acts to 
expand a dense cluster of drops into the hot ambient gas, thus promoting 
evaporation. For dilute clusters of drops, neither electrostatic nor turbulence 
effects were important. Based upon single-droplet-stream experiments by 
Sangiovanni and Liscinsky 16 and the results of the calculations, inferences were 
made about how electrostatic spray dispersion may reduce the soot emission index 
by increasing drop spacing. Also, limited results pertaining to the ignition of 
nearly-dense clusters of electrostatically charged drops were discussed. 

The case of a jet flow of charged drops has been formulated, but has not 
been programmed for numerical calculations. This would require a subsequent 
study to obtain information on dispersion and evaporation effects in a geometry 
more 1 i in our experiments. 
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ELECTROSTATIC DISPERSION OF DIESEL FUtu JETS AT HIGH BACK PRESSURE 

+ ++ + 

E. Y. Kwack, L. H. Back and C. P. Bankston 


Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, CA 91109 


ABSTRACT 


An experimental study of electrostatically atomized and dispersed fuel jets 
has been conducted for various back pressures to 600 psig. No. 2 diesel fuel was 
injected through an electrostatic spray triode designed for a small diesel 
engine. Charge density measurements were conducted at various combinations of 
injection velocities, electric potentials and back pressures. The charge density 
of fuel drops linearly increased up to 1.5 C/m 3 with increasing electric 
potential until breakdown occurred. After breakdown the charge density was 
reduced by 40 to 60 % and again increased with a lower slope as electric 
potential increased. At higher flow rates, breakdown occurred at higher 

voltages. At higher back pressure, lower charge density was obtained and 
breakdown occurred at higher voltages. Visual observation showed that 
significant electrostatic dispersion was accomplished at high back pressures, and 
that the average drop size was about the same as the spray triode orifice 
diameter. 


I. INTRODUCTION 

Diesel engines are widely used for transportation vehicles and power 
generation. Unfortunately, one of their important drawbacks is the fact that 
they produce non-negl igible amounts of soot. The soot is formed in the fuel rich 
regions of the spray where a large amount of fuel vapor accumulates without 
burning and undergoes pyrolysis reactions which transforms it to soot. Mixing 
and controlled atomization have good potential to reduce soot yield in spray 
flames (Prado et al, 1977). Many methods for improved mixing and dispersion of 
the spray to reduce local fuel rich regions have been investigated for improved 
combustion and reduced soot formation. 

In general, the breakup of fuel jets in diesel engines is accomplished by 
high pressure injection in the range 20 MPa, 200 atm (for indirect injection) up 
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to about 140 MPa, 1400 atm (for direct injection). The mechanism of break up of 
a fuel jet at high velocity injection has been investigated and was well reviewed 
in articles by Reitz and Bracco (1982) and Lin and Kang (1987). Swiri is 
imparted to the flow field in the chamber to further assist the fuel -air mixing 
process, and soot formation is strongly dependent upon swirl characteristics. 
However, fuel rich regions apparently still persist in current engines, and 
methods for improved mixing and reduced soot formation continue to . be of 
interest. 

Techniques for improving in-cylinder combustion characteristics include 
mechanical (Pischinger and Cartellieri, 1972; Khan et al , 1972), chemical 

(Gaffney et al , 1980) and electrical methods. Electrical concepts involve the 
injection of a fuel jet which is electrostatically broken up and dispersion 
enhanced by the repulsion of like-charged droplets. Numerous attempts have been 
made to develop electrostatic fuel atomizers for use in combustion systems and 
also to develop models to explain the electrostatic processes of breakup and 
dispersion of fuel jets, and those investigations were reviewed by Kelly (1982). 
However, electrostatic dispersion of a fuel spray has not proven practical due to 
the very high resistivities of liquid hydrocarbon fuel (10^ - 10^ ohm-m) 
(Kelly, 1981, 1982; Luther, 1962) which presents difficulties in producing 

electrostatically charged sprays at flow rates above 0.1 m£/s (Luther, 1962; 
Newgard and Noon, 1964). The electrostatic spray device conceived and 
demonstrated by Kelly (1981, 1982) has overcome previous limitations and can 
charge mineral oil Marcol-87 (viscosity * 25cp) up to 1.8 C/m^ for flow rates up 
to 1.25 ml/s. The detailed schematic diagram of the spray triode can be found in 
Refs. (Kelly, 1982; Bankston et al. 1988). Bankston et al (1988) used a spray 
triode provided by Kelly and obtained charge densities of 1.5 - 2.0 C/m^ for a 
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mineral oil and diesel fuel oil. Of note is that both the experimental studies 
by Kelly (1982) and Bankston et al . (1988) were conducted at ambient back 

pressure. No experimental study has been conducted with this kind of triode for 
high back pressures. 

Bankston et al . (1988) also estimated the effect of charged droplets on the 
dispersion of a spherical cluster of fuel drops under diesel engine conditions 
using a model proposed by Bel lan (1983) based on a model of evaporation of dense 
clusters of drops (Bel lan and Cuffel, 1983). That model predicts the expansion 
of an initially spherical droplet spray cloud of radius R 0 , due to the repulsion 
of uniformly charged and distributed droplets. Bankston et al (1988) has 
calculated spray cloud dispersion versus charge density in 0.5 ms; for example, a 
spherical spray volume 1 cm in radius and containing 54 mg of diesel fuel will 
expand by a radius factor of 1.4 in 0.5 ms at a charge density of 2 C/m 3 . 
Although the model does not account for drag effects and is developed for a 
spherical spray cluster, not for a cylinderical plume like a spray jet from ? 
nozzle, the results indicate that the electrostatic dispersion can be considered 
promising to control or reduce soot formation. A simple model of electrostatic 
spray plumes was developed by True (1983) but no account was taken of droplet 
evaporation and droplet interactions due to their near proximity. Harstad and 
Bel lan (1988) proposed a more complex model for a spherical cluster of charged 
drops by extending their model for evaporation of dense Glisters of fuel drops 
(Bell an and Harstad, 1987, 1988) to account for electrostatic effects. 

Electrostatic spray concepts have been applied to the actual combustion 
chamber (Ahmad et al, 1980; Miller et al, 1987) however, there was a lack of 
study of the fundamental phenomena. Moreover, not much experimental work has 
been done at high back pressure. With this as a background, this detailed 
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experimental study was carried out to investigate the effects of high back 
pressure on the charging of fuel jets and the behavior of breakup and dispersion 
of charged fuel drops. 

II. EXPERIMENTAL SF.TUP 

The nozzle used in the current study is the first electrostatic spray triode 
(EST) "single hole" injector designed by Parker Hannifin Company to fit into a 
Volkswagen diesel four-cylinder engine of nominal 48 horsepower output. The 
basic operational principle of the EST nozzle can be found in previous reports 
(Kelly, 1982; Bankston et al., 1988). The OMC (Oxide-Metal Composites) company 
manufactured the insulator-electrode (emitter) assembly and Parker Hannifin 
Company assembled it into the housing. An;, diesel fuel or fuel simulant can be 
used in this EST nozzle as long as the fluid is non-conducting and has similar 
electrical properties as diesel fuel. The EST nozzle has a 254 pm diameter 
orifice and can be operated up to fuel pressures of 3000 psig, electric 

potertials of 20 kV and temperatures of 200°C. 

The EST was located in the top of a high pressure chamber 15 cm I.D. and 40 
cm height which also has a 10 cm diameter window to allow spray development to be 
observed (Fig. 1). The EST injects fuel vertically downward into a cylinderical 
glass container of 10 cm diameter and 25 -cm height placed inside the high 
pressure chamber. The glass container is lined with copper wire screen on the 
wall and has copper gauze on the bottom in order to collect the charge of the 
dispersed fuel spray. 

The copper wire screen and gauze are connected to an electrometer to measure 
the electric charge carried by the fuel per unit time if which when divided by 
the volume flow rate 0 gives the mean charge density p e . 
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Da = if/Q 


( 1 ) 




The anode is connected to the electrometer as shown in Fig. 1 by a switching 
arrangement when the total current is measured. 

The high pressure chamber was pressurized up to 600 psig by N£ gas. Fuel is 
supplied to the EST nozzle by the pressure difference between the fuel tank and 
the high pressure chamber through a line with a 10 pm filter to remove particles 
which can cause the nozzle to plug, to spray erratically, and to lead to early 
electrical breakdown. The flow rates through the EST were calibrated as a 
function of pressure difference measured using a differential pressure transducer 
for three different chamber pressures of 0, 300 and 600 psig. Ic was found that 
the calibration was almost independent of the chamber pressure. Since it was 
planned to increase the pressure ; nside the chamber to 600 psig to simulate the 
pressure inside a small diesel engine cylinder before the injection of fuel, the 
high pressure chamber and the fuel tank were located in one room for safety and 
the associated instrumentation was located in an adjacent control room. For 
remote control, a solenoid valve was located in the fuel line between the 10 
filter ana the EST nozzle. 
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The test fluid was No. 2 diesel fuel with a viscosity of about 6cp. The 
density and the surface tension of No. 2 diesel fuel are 0.824 g/cm 3 and 30 
dyne/cm, respectively. The electric charge density of the diesel spray was 
measured for various flow rates/injection velocities at various chamber pressures 
of 0, 300 and 600 psig by applying different electric potentials with a power 
supply. 

Characterization of the spray was carried out by flash illumination and 
photographs. Analysis of the photographs was carried out using an image 
processing system. A photograph was placed on a light table in view of a video 
camera (Model CC-65S, DAGE-MTI, Inc). The analog signal from the video camera 
was fed to a digitizer which writes the data into RAM memory using an add-on 
board (PCvision Frame Grabber, Imaging Technology) in conjunction with an IBM PC- 
AT. A computer code was developed with the aid of commercially available 
software ( IMAGELAB/IMAGET00L, Werner Frei Associates) for image enhancement and 
analysis of droplet size. 

III. RESULTS AND DISCUSSION 

1 . Electric Charge Density 

The amounts of charge emitted by the cathode and that swept downstream per 
unit time, if and if, respectively, were measured for a wide range of the flow 
rate at p c « 0 psig by a switching arrangement. The results of if and if for a 
fixed flow rate Q - 1.2 mi/s are shown in Fig. 2 as a function of the electric 
potential. The cathode starts emitting electrons from its surface if the local 
electric field is higher than the minimum value needed to extract electrons from 
metallic surface which is of the order of 10® V/m. The number of electrons 
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emitted from the cathode per unit time is proportional to the strength of the 
1 al electric field above the minimum value at a given fuel flow rate (Fig. 2). 

Once the electrons are injected into the fuel, they are coupled with the 
fuel near the cathode, forming a charged cloud or layer. Some of the charged 
fuel is attracted to the anode/orifice electrode (i a ) and the other is discharged 
to the collector (if) as shown in Fig. 2. The amount of charge attracted to the 
orifice electrode depends on electron mobilities of the charged fuel which are in 
the range of 10' 7 to 10‘ 8 m 2 /V.sec (Inuishi, 1979). Therefore, more charged fuel 
is attracted to the orifice electrode as the electric field increases (Fig. 2). 
If the fuel velocity or flow rate increases at a fixed electric potential, more 
charged fuel in the layer moves downstream (Fig. 3) and the charged layer becomes 
thinner which results in more emission of electrons from the cathode. In other 
words, the total charge emitted from the cathode increases with increasing flow 
rate (not shown) like the charge carried downstream does (Fig. 3). 

The layer or the charged fuel becomes thicker as the electric field 
increases, until it reaches the anode locally and forms an electric channel or 
low conducting path, so-called electric breakdown. After elecric breakdown 
occurs, the amount of electrons emitted by the cathode increases drastically by 5 
to 10 times (Fig. 2). However, since most electrons move to the anode directly, 
and are not carried by the fuel, the amount of electrons emitted by the cathode 
no longer depends on the fuel flow rate. At the same time, most of the electrons 
are attracted to the anode and only, a small amount is swept away by the fuel so 
that if is reduced by 40 to 60 % after breakdown (Figs. 2 and 3). The ratio of 
the charge swept downstream to the total emitted charge was observed to be 0.1 to 
0.5 for the range of flow rates 0.4 to 1.8 mi/s before breakdown occurred. After 
breakdown, the ratio was reduced below 0.05. After breakdown if and if again 
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increase gradually as electric potential increases. However, the slope of if is 
much lower than that before breakdown while the slope of if remains the same. 
Therefore, the ratio decreases as electric potential increases after breakdown. 
Since the rate of growth of the charged layer decreases as the fuel velocity 
increases, breakdown occurs at higher voltages, -8.5 kV at 0.51 mi/s to -13 kV at 
2.5 mi/s (Fig. 3). 

The maximum amount of charge on a liquid drop is limited by the Rayleigh 
relation: 


p e * 12(2cq7) 1 / 2 d" 1 • 5 (3) 

max 

If the charge increases more than the maximum limit, the liquid drop would 

breakup into several small drops. Therefore, the charge density p e is an 

max 

important parameter in determining the largest drop size. 

The experimental results shown in Fig. 3 were converted to charge density by 
dividing if by the flow rate Q and the results are shown in Fig. 4. The charge 
per unit time carried by tho fuel to the collector, if, increased with increasing 
flow rate to 2 mi/s at a fixed voltage before breakdown occurred (Fig. 3). 
However, the incremental increase was not enough to give the same or larger 
charge density, so that the charge density monotonically decreased with 
increasing flow rate. However, the maximum charge density which was obtained 
just before electrical breakdown occurred is 1.0 C/m 3 at the flow rate of 0.51 
mi/s and increases with increasing flow rate, reaching a peak value of 1.53 C/m 3 
at a flow rate of 1.5 mi/s. Further increase of the flow rate resulted in a 
slight decrease of the maximum charge density. 

Since the chamber pressure in a diesel engine before fuel injection is 
relatively large, the charge density of No. 2 diesel fuel was measured for high 
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back pressures up to 600 psig (at 300 and 600 psig) and the results are shown in 
Figs. 5 and 6. The results for the high back (chamber) pressures exhibit similar 
behavior as those at p c * 0 psig i.e. the charge density increases with an 
increase of the applied voltage at a fixed flow rate, and decreases with 
increasing flow rate at a fixed applied voltage. The chamber pressure, however, 
plays an important role in determining the charge density and the breakdown 
voltage. Breakdown seems to occur at a higher voltage for a fixed flow rate as 
the chamber pressure increases. For example, for Q=0.51 mi/s the breakdown 
occurs at E=-8.5 kV for p c = 0 psig, at E *-10 kV for p c = 300 psig and at E=-ll 
kV for p c = 600 psig. The decrease in the charge density right after breakdown 
is much smaller for higher chamber pressures. The charge density decreases by 
20% for p c * 300 and 600 psig, compared to 50% for p c * 0 psig at Q = 0.51 mi/s. 
Moreover, the charge density decreases as the chamber pressure increases before 
breakdown occurs, if the flow rate and the applied voltage are the same. This 
behavior is clearly seen in Fig. 7 which shows the charge density as functions of 
the flow rate at E = -8 kV for p c * 0, 300 and 600 psig. The charge density at 
p c * 0 psig is about 1 C/m 3 at Q * 0.51 mi/s and monotonically decreases with 
increase of the flow rate, reaching about 0.3 C/m 3 at Q = 2.5 mi/s. The charge 
densities at pressures of 300 and 600 psig show similar behavior but, however, 
are smaller than those at 0 psig by 10 - 17% for 300 psig and 40 - 50 % for 600 
psig, depending on the flow rate. 

The maximum charge densities measured at three different back pressures are 
plotted in Fig. 8 using the results shown in Figs. 4, 5 and 6. This plot clearly 
demonstrates that the maximum charge density becomes smaller as the back pressure 
increases if the electric potential is limited, i.e. E - -14 kV. The maximum 
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charge density is smaller by 20 - 26% at 300 psig and by 22 - 43% at 600 psig 
. _ rnmnarfid to that at 0 psig* 

depending t?fi i iGw » wL, " r nf 

The detailed maximum charge density with the high pressure spray triode was 
measured at two different chamber pressures p c of 0 and 300 psig for a wider 
range of flow rates of 0.3 to 4.2 mi/s. In this case, the electric potentiai was 
applied to -18HV. The maximum charge density increases with increasing the flow 
rate, reaching maximum values of 1.65 C/m» at p c • 0 psig and 1.5 C/m» at p c - 
300 psig as shown in Fig. 9. It remained at the maximum value for a range of 
flow rates and tended to decrease with further increase of flow rate. The 
maximum charge density at 300 psig was lower than that at 0 psig by 10-25%, 

depending on flow rate. 

The electrostatic spray triode used in the current study shows different 
characteristics at low flow rates compared to those used by Kelly (1982) and 
Bankston et al (1988) for atmospheric pressure. For example, Kell, reported 
almost constant charge density of 1.5 C/m» for a range of flow rates from 0.25 to 
1 25 mi/sec. However, the charge density in Fig. 9 almost increases linearly up 
to flow rates of 1.5 mi/s at 0 psig and 2.2 mi/s at 300 psig, and then remains 


constant. In any case, the data confirm that charge densities of 1.5 C/m* up to 
a flow rate of 4.3 mi/s are definitely achievable at atmospheric pressure, and 
slightly lower values are found at high pressures. The performance at g 
pressure may be improved by redesigning the insulator which limits 
potential. If an electrostatic spray triode were used in a diesel engine, the 
applied voltage should be varied according to the injected fuel amount or 
injection velvety since the maximum charge density is obtained at higher 
electric field as the flow rate increases. 
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2. Flow Visualization 


Pictures were taken of diesel fuel injected through the high pressure spray 


t v»*l rirln 
I I uuv, « 


p^r>+os were taken through the 10 cm I.D. window at various combinations 


of test conditions of chamber pressure, fuel injection velocity and applied 
voltages as shown in Table 1. 



Table 1. 


pressure, psig 

V i * 20 m/s 

V i = 40 m/s 

0 

0, -8, -lOkV 

0, -8, -lOkV 

300 

0, -10, -12kV 

0, -10, - 12k\ 

600 

OkV 

OkV 


Without Electric Charge 

When fuel was injected into air at 0 psig (atmospheric pressure) a single 
stream of fuel was observed all the way across the window without electric charge 
even at an injection velocity of 40 m/s (Fig. 10a). At p c = 300 psig, = 20 
m/s, the fuel jet breaks up even without the applied voltage since the nitrogen 
density is about 21 times larger than that of air at atmospheric pressure (Fig. 
10b). The fuel jet begins breaking up at 0.75 cm beyond the injection nozzle 
(xj/d n * 30) and mixes with surrounding nitrogen gas, f orming a spray plume 
having an included angle of 26°. However, a dense central core exists for 
several centimeters beyond the injection nozzle. Many large irregular droplets 
were observed to remain in this dense core region. When fuel is injected at 40 
m/s (Fig. 10c), the jet breaks up earlier, less than 0.5 cm beyond the nozzle 
exit, forms a narrow spray angle of 23°, and breaks up into much smaller 
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droplets. The fuel was also injected into nitrogen gas having a pressure of 600 
psig to simulate an actual diesel engine. Compared to data at 300 psig, spray 
angles are larger and the jet. breaks up earlier: 6 - 35°, x x - 0.5 cm at 20 m/s 
(Fig. lOd) and 6 - 30°, = 0 cm at 40 m/s (Fig. lOe). The droplet sizes are 
somewhat larger at 600 psig if the injection velocity is the same. Despite the 


large spray angle at p c - 600 psig, the dense core still exists for several 

centimeters beyond the injector nozzle. 

Atomization of fuel jets has been studied using aerodynamic interaction 


mechanisms (Ranz, 1958; Reitz and Bracco, 1982; Wu et al, 1986). The following 
relationships were proposed for initial spray angle a, and initial average drop 


diameter d: 


tanC^/ 2 ) “ 1- 4* ( ^ J 1 / 2 

A Pf 

d • B x max( T ) 

PgVi 2 

where the aerodynamic breakup parameter T • pf/pg (Re/Wef) 2 . The functions f and 
x max tend to 3 l / 2 /6 and 3/2, respectively when T is greater than unity. A and B 
are constants whose values depend on the nozzle geometry and must be determined 
experimentally. Of note is that the Weber number Wef is based on the fuel 
properties, injection velocity and on nozzle diameter, not on the qas properties 
and on fuel droplet size. 

The spray angle is essentially linearly proportional to the square root of 

gas density, being independent of injection velocity if the function f is 

constant. This relation for spray angle was confirmed by extensive experiments 
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(4) 

(5) 


I 


; by Reitz and Bracco(1982) . In the current experimental results, the Reynolds 

numbers are 700 at 20 m/s and i4uu at 40 m/s. I he values of the aerodynamic 
breakup parameter T are 2 at 20 m/s, a .^d 0.5 at 40 m/s when p c = 300 psig. Since 
the density of nitrogen at 600 psig is 1.95 times larger than at 300 psig, those 
at 600 psig are about one half of those at 300 psig, being 1 and 0.25, 
respectively. Ihe values of the function f from Reitz and Bracco(1982) were 
calculated to be 0.21 at V, - 20 m/s and 0.18 at V 1 - 40 m/s when p c - 300 psig; 
and 0.20 at Vi - 20 m/s, and 0.18 at Vi - 40 m/s when p c - 600 psig. The values 
of the geometrical parameter A for the electrostatic spray triode were calculated 
to be 2.02, 1.99, 1.96 and 2.05 with an average of 2. The current data support 
that the spray angle of the fuel jet can be predicted if the geometrical 
parameter A is known for a particular nozzle. 

The average droplet size from Eq. (5) is inversely proportional to the 
square of the injection velocity and to the density of gas which the fuel is 
injected into, if x max is constant. Although the detailed sizes of droplets were 
i not measured at E » 0 kV, Fig. 10 clearly shows that the droplet size becomes 

much smaller if the injection velocity is doubled at the same pressure. However, 
when the injection velocity is the same, the droplets tend to increase in size 
| with increase in gas density or pressure which is opposite to Eq. (5). A recent 

experimental study by Wu et al. (1986) reported similar behavior. Drop 

| coalescence leads to a rapid increase in drop size, particularly at high gas 

density, which may be prevented by electrostatic repulsion. 

i 

i With Electric Charge 

When fuel was injected into air at 0 psig (atmospheric pressure) a single 
j stream of fuel was observed all the way across the window without electric charge 

I 

i 


even at the flow injection velocity of 40 m/s (Fig. 10a). However, the fuel jet 
i m uhon tho unitanp i c annlipH. not shown since similar photos were shown 

ISI WWirw U fs nnvit vtiw ■ x' • * — *T T * * / 

in an earlier paper (Bankston et al , 1988). As the charge density increased, the 
break up of fluid occurred earlier, and the spray angle increased more rapidly. 
It seems that the droplet size becomes smaller as the charge density increases. 

At p c - 300 psig, and Vi » 20 m/s the fuel jet breaks up without the applied 
voltage (Fig. 10b). However, the applied voltage caused the fuel stream to break 
up considerably more and increases the initial spray angle to almost 180° at E ■ 
-12kV (Fig. 11a). A large number of small droplets were observed everywhere 
except in the dense core which extended to about 1.5 cm beyond the nozzle at E -• 
-10kV(Fig. lib). At E - -12kV, ?hr>re was significant fuel jet breakup in the 
near proximity of the orifice, and large lateral dispersion of the charged 
droplets. If soot is primarily formed in the fuel rich region such as a dense 
core region, electrostatic dispersion should reduce the formation of soot. 
Similar behavior was observed for the higher injection velocity - 40 m/s. 
Since the inertia force is larger, however, the dense core region extended to 3 
cm at E - -lOkV (Fig. 11c) and to 2.5 cm at E - -12kV (Fig. lid). 

The drop size distributions were obtained using image processing for Figs. 
11a, lib and lid. Details on the procedure are given in the Appendix. The drops 
in the three areas in Fig. lib were examined. About 20 drops in each region were 
processed and the average diameters, djo* were estimated to be 235, 235 and 240 
(i, respectively. Since the drop size distributions in the three regions are 
similar, the values of djQ and ^32 for all drops in the three regions were 
calculated. They are 236 and 242 **, respectively. The detailed distribution is 
shown in the bottom of Fig. 12. The average charge density of fuel drops was 
obtained from Fig. 5 and was plotted in Fig. 13 as a function of the drop 
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are included as solid lines. The current averaged value is higher than Kelly's 
prediction but is lower than the Rayleigh limit. This is the same as the 
previously reported result (Bankston et al . , 1988). The same estimation of drop 
size was done for one region in Fig. 11a. The drop size distribution is shown in 
the top of Fig. 12. Even though the charge density is higher than the previous 
case (Fig. 5), the drop sizes are larger. The values of dig and d 32 are 367 and 
384 n , respectively. As shown in Fig. 13, the data is close to the Rayleigh 

1 imit. 

In the case of the higher injection of V-j = 40 m/s (Fig. lid) the drop size 
distribution is similar to that of = 20 m/s at E = - 10 kV even though the 
aerodynamic force on the fuel jet is four times larger. The values of d 10 and 

d 32 are 234 and 241 n, respectively (sample size, Nj 0 t a l = 85 )- The char 9 e 

‘ density was about 0.9 C/m 3 which is the same order of that of V, = 20 m/s at E = 

j -10 kV ( pq = 0.78 C/m 3 ). Therefore, it seems that the drop size is determined by 

the charge density rather than by the aerodynamic forces. Since high injection 
J of fuel tends to increase the fuel rich or dense core region, better performance 

is expected at the lower injection velocities. However, the performan of 
I electrostatic dispersion should be optimized by adjusting various parameters such 

| as the injection velocity, electric potential and nozzle size to inject a certain 

amount of fuel in an engine to produce a suitable size of drops and less fuel 
j rich or dense core regions. A flow visualization study was attempted for the 

higher back pressure of 600 psig, but it was not successful due to the failure of 
I the insulation in the electrostatic spray triode. 

I 

I 

I 
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Summary and Conclusions 

The results of this investigation are summarized as follows: 

1. Charge density measurements (p c = 0 psig) 

i. The total charge emitted by the cathode per unit time increased 
linearly with the electric potential until breakdown occurred. After 
breakdown occurred, the total current sharply increased by 5 to 10 
times depending on flow rates and again gradually increased as the 
electric potential increased. 

ii. The charge swept downstream by the fuel showed a similar behavior as 
the total charge, but decreased by 40 to 60 % depending on flow rate 
when breakdown occurred. 

iii • The ratio of the charge swept downstream to the total emitted charge 
was 0.1 to 0.4 before breakdown and was smaller than 0.05 after- 
breakdown. 

iv. As the flow rate increased, breakdown occurred at higher voltages; -8.5 
kV at 0.51 mi/s to -13 kV at 2.5 mi/s. 

v. The maximum charge density increased with flow rate , reached a maximum 
value of 1.65 C/m-\ remained at the maximum value for a range of flow 
rates and decreased slightly with further increase of flow rate. 

2. Charge density measurements (p c « 300 and 600 psig) 

i. The charge density at high back pressures exhibited similar behavior as 
those at p c - 0 psig; the charge density increased with electric 
potential at a fixed flow rate and decreased with flow rate at a fixed 
applied voltage. 
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ii. The charge density decreased with increasing back pressure before 
breakdown occurred; being smaller than that at 0 psig by 10-17 % for 
300 psig and by 40-50 % for 600 psig when Q = 0.51 mZ/%. 

iii . The maximum charge density was also smaller by 20-26 % for 300 psig and 
by 22-43 % for 600 psig depending on flow rate when the applied voltage 
was limited below -14 kV. 

iv. Breakdown occurred at a higher voltage as the back pressure increased; 
at -8.5 kV for p c = 0 psig, at -10 kV for p c = 300 psig and at -11 kV 
for p c = 600 psig when Q = 0.51 mi/s. 

v. The decrease in the charge density right after breakdown was much 
smaller at higher back pressures; 20 % reduction for 300 and 600 psig, 
compared to 50 % for 0 psig at Q = 0.51 mi/s. 

3. Flow Visualization 

i. In general, with electric charge, rapid breakup of the jet was achieved 
without large droplets in the core region. 

ii. The spray angle increased to larger than 180° with increase of 
electric potental which enhanced dispersion of fuel near the nozzle 
exit and reduced the dense core region. 

iii. It was observed that the size of drops may be determined by the charge 
density rather than by aerodynamic forces. 

The present investigation has indicated significant early breakup and 
dispersion of charged diesel fuel jets by using an electrostatic spray triode at 
high back pressures compared to aerodynamic breakup and dispersion. These 
results are encouraging for engine injection, but demonstration of this technique 
under pulsed fuel injection and autoignition conditions requires subsequent 
investigations to evaluate the potential for improvement in combustion and 
reduced soot formation. 
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Nomencl ature 

A jeometrical constant in Eq. (5) 

B geometrical constant in Eq. (6) 

d drop diameter, m 

d n spray nozzle diameter, m 

djo arithmetic mean diameter, m 
d 32 volume- surface mean diameter, m 

E electric potential across cathode and anode electrodes, V 
f function in Eq. (4) 

i a charge attracted to anode per unit tin; amp 
if charge swept downstream by fuel per unit time, amp 

it total charge emitted by cathode per unit :inu-, amp 

mf fuel mass, g 

N number of fuel drops 

p c chamber pressure, psig 

Q fuel flow rate, mi/s 

Re Reynolds number, pfV-j d n /p^ 

T aerodynamic breakup parameter 

fuel injection velocity, m/s 
Wef Weber number based on fuel properties, pf^V-j 2 d n/7 

x max function in Eq. (5) 

xj length of intact core, m 

7 surface tension of fuel, N/m 

c 0 permittivity of free space, 8.85 x 10 -12 C/V-m 
e spray angle, degree 
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/if viscosity of fuel, kg/m-s 
Pe charge density, C/m 3 
°e,max maximum charge density, C/m 3 
pf fuel density, g/cm 3 
pg gas density, g/cm 3 
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Appendix 


Image Processing 

Drop sizes were determined in the following manner. Image frames with 512 x 
480 pixel resolution were formed over areas approximately 2 cm on a side from a 
spray photograph. The corresponding pixel size (image physical resolution) is 
about 45 /im. Any drops of about this size or smaller would not register in the 
image; at most, they would add an apparent "noise" to the image. Thus 
measurements were valid down to drop sizes of order 100 /im. Using a mouse guided 
cursor, isolated, distinct drops in a frame were selected. Pixel values along 
two perpendicular lines through a drop center were obtained. The first and 
second moments of these pixel arrays were used to determine an average pixel 
threshold value corresponding to the RMS drop extension in each direction. The 
image in a box of pixels surrounding the drop was thresholded and a pixel count 
of the drop in the threshold image made. Assuming a Gaussian pixel gray level 
distribution for a drop, the actual drop cross-sectional area is calculated as 
twice this pixel count. 


of Figures 

Schematic diagram of experimental setup 

Currents vs electric potential at p c « 0 psig and Q = 1.2 mi/s 

Charge swept downstream by fuel per unit time vs electric potential for 
various flow rates at p c 0 psig 

Charge density vs electric potential for various flow rates at p c - 0 psig 

Charge density vs electric potential for various flow rates at p c * 300 psig 

Charge density vs electric potential for various flow rates at p c - 600 psig 

Charge density vs flow rate for various back pressures at E ■ -8 kV 

Maximum charge density vs flow rate at back pressures of 0, 300 and 600 psig 
and at E , -14 kV 

Maximum charge density vs flow rate at back pressures of 0 and 300 psig and 
at E = -18 kV 

Diesel fuel jet spray at E = 0 kV, = 40 m/s and p c = 0 psig 
Diesel fuel jet spray at E = 0 kV, * 20 m/s and p c = 300 psig 

Diesel fuel jet spray at E - 0 kV, » 40 m/s and p c = 300 psig 

Diesel fuel jet spray at E = 0 kV, = 20 m/s and p c = 600 psig 

Diesel fuel jet spray at E = 0 kV, V, = 40 m/s and p c = 600 psig 

Diesel fuel jet spray at E = -12 kV, Vj = 20 m/s and p c = 300 psig 

Diesel fuel jet spray at t = -10 kV, Vj = 20 m/s and p c = 300 psig 

Diesel fuel jet spray at E - -10 kV, Vj = 40 m/s and p c = 300 psig 

Diesel fuel jet spray at E ■ -12 kV, V, = 40 m/s and p c = 300 psig 

Droplet size distributions of fuel sprays of Figs. 11a and lib 
Charge density 's droplet size 
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Abstract 


A theory of evaporation and dispersion of electrostatically charged drops 
has been developed for drops belonging to a spherical cluster exposed to a flow. 
Under the assumption of constant atmospheric pressure, the quasi -steady 
approximation was made for the gas phase whereas the drop-temperature history is 
unsteady. The model takes into account interdrop interactions (in terms of heat 
and mass transfer) due to drop proximity, turbulence exchange processes between 
the cluster and its surroundings and electrostatic force effects due to the 
charge on the drops. Calculations based upon this model were made fer charged 
as well as uncharged clusters of drops. The charge was varied from a null value 
to the maximum possible charge found empirically for hydrocarbon sprays. 
Moreover, the turbulence model was varied in such a way as to simulate the 
cluster embedded into a flow where turbulence develops with time (Model 1) or a 
flow with pre-existing turbulence (Model 2). The results show that the control 
parameters for the evaporation of charged drops are different from those for 
uncharged drops in dense clusters; turbulence levels which were shown to be 
crucial for the latter in the dense cluster regime do not affect the former in 
the same regime. For dilute clusters turbulence is unimportant in both cases. 
Moreover, drop charging does not affect dilute clusters of drops whereas dense 
clusters of drops are substantially affected. Based upon existing experimental 
data, inferences are made about how electrostatic spray dispersion can affect 
soot control in power systems using fuel sprays. Limited results pertaining to 
the ignition of nearly-dense clusters of electrostatically charged drops are 
discussed as well . 


I. Introduction 


Power systems using fuel sprays are widely used for mobile and stationary 
power generation. Unfortunately, one of their important drawbacks is the fact 
that they produce non-neg' : aible amounts of soot. The soot is formed in the 
fuel rich regions of the spray through the following process. In general the 
fuel rich regions of the spray correspond to dense clusters of drops evaporating 
and producing a large amount of fuel vapor that is not consumed fast enough in 
the flame. This fuel vapor accumulates without burning and thus instead of 
being transformed into the usual products of combustion undergoes another set of 
reactions (pyrolysis reactions) which transform it into soot. Consistent with 
this observation, mixing and controlled atomization have good potential to 
reduce soot yield in spray flames (Prado et al., 1977). 

Among various means of controlling soot formation, electrostatic charging 
of the drops was considered promising by Bellan (1983) in the particular case of 
Diesel engines. Previous investigators (Mayo and Weinberg, 1970; Thong and 
Weinberg, 1971) had also considered dispersing particles using electric effects, 
although not specifically for Diesel engines. Mayo and Weinberg (1970) have 
shown that an applied electric field can reduce both the size and the number 
rate of formation of soot particles in atmospheric diffusion flames. The 
situation in Diesel engines was also investigated by Ahmad et al. (1980) and it 
was shown that, unlike in Mayo and Weinberg (1970), the mass of particulates 
emitted increased substantially when electric fields were applied. The 
difference between the results of Mayo and Weinberg (1970) and Ahmad et al . 
(1980) was attributed to the basic difference of the two combustion systems, but 
no physical mechanisms accounting for the observations could be provided because 
there was a lack of study of the fundamental phenomena. In Bellan (1983), 
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through a very simple calculation, based upon a model of evaporation of dense 
clusters of drops (Bel lan and Cuffel, 1983), it was shown that electrostatic 
drop charging could be achieved in a Oiesel engine and that it would disperse 
the drops from a dense configuration to dilute configurations where soot 
formation is no longer a problem. Moreover, the rdded benefits of secondary 
atomization induced by charging were also mentioned in Bellan (1983). Roth and 
Kelly (1983) showed indeed that convulsive disruption occurs during charged- 
droplet evaporation, yielding a limited number of siblings, about seven. 
However in Bellan (1983), no definite conclusions could be asserted because the 
model used was rather restrictive in its assumptions. A simple model of 
electrostatic spray plumes was also developed by True (1983), but in that model, 
unlike in Bellan (1983), no account was taken of droplet evaporation or droplet 
interactions due to their proximity. However, limited drag effects were 
accounted for by True (1983) unlike the model upon which the work of Bellan 
(1983) was based ( Bellan and Cuffel 1983). 

Recent progress in modeling of dense collections of evaporating drops 
(Bellan and Harstad, 1987a, 1987b) and also new experimental and theoretical 
results related to electrostatic atcmization (Sangiovanni and Liscinsky, 1984) 
have reopened the question of whether electrostatic spray atomization is a good 
way to disperse drops and if so in which regimes should it be considered. 
Additionally, one should know what are its disadvantages, if any. 

The present work answers some of these questions and shows that the 
control parameters for the evaporation of electrostatically charged dro^s 
clustered together are different from those of uncharged drops in clusters in 
the dense cluster regime. The present model pertains to combustion at 
atmospheric pressure. Instead of performing soot calculations which could be 
limited by the use of a simple kinetic model, we make Inferences regarding the 
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impact of charging upon soot control in heterogeneous flames by using already 
existing experimental observations of soot formation in well-defined sprays 
(Sangiovanni and Liscinsky, 1984). Since this model is expected to yield 
predictions of qualitative rather than quantitative value, no attempt is made to 
make detailed comparisons of experiments with the idealized situation considered 
below. I he focus is here on showing that electrostatic spray charging can 
conceivably be used to disperse drops in a regime that is useful to control soot 
production. 

II. Formulation 

Under consideration is a monodisperse, uniformly distributed cluster of 
drops in a spherical configuration where each drop may be equally charged. The 
cluster as a whole is assumed to move in a straight line trajectory through an 
ambient gas at constant pressure, and as it does so, the cluster may undergo 
radial expansion due to mutual repulsion of the charged drops. 

The gas inside the cluster is at the same, constant, pressure as the gas 
outside the cluster. On a large length scale (many drop diameters), the cluster 
is taken to be spatially homogeneous in thermodynamic quantities. A coordinate 
system fixed with the ambient gas is used. In this system, the drops hav' two 
velocity components as shown in Figure 1: a uniform axial component along the 
trajectory direction and a radial expansion component which is assumed to be 
self-similar. The gas inside the cluster has an Initial null velocity (like the 
ambient gas) but It evolves to have axial and radial components In the same 
manner as the drops. For a very dense cluster, momentum exchange between the 
drops and gas Inside the cluster will rapidly reduce the relative (slip) axial 
velocity between the two (Bellan and Harstad, 1987a). The cluster will then act 
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as one large solid entity in its movement in the ambient. Note that this does 
not mean that radial slip is negligible. Calculation of mass loss for uncharged 
clusters or ambient gas ingestion for charged clusters, which are dependent on 
radial slip, is an important aspect of the model. 

Previous results (Bellan and Harstad, 1988) have pointed out that turbulent 
exchange between the cluster and the surroundings are crucial during evaporation 
of dense clusters of drops. In the present model turbulent exchange processes 
are considered as well and their relative importance with respect to 
electrostatic effects is of interest. 

Similar to previous models (Bellan and Cuffel, 1983; Bellan and Harstad, 
1987a and 1988), here each drop of the cluster is considered surrounded by a a 
fictitious sphere of influence centered at the drop's center and having a radius 
which is the half distance between two adjacent drop centers. Thus the model 
consists of two coupled submodels: that of an individual drop in its sphere of 

influence and that of the multitude of interacting spheres of influence inside 
the cluster, as follows. 

A. Model q £ one drop inside its sp here of influpnrp 

The classical approach is taken here whereby the radial variation of gas 
properties in a sphere of influence is analytically estimated using quasi-steady 
equations without convection, and the enhancement of evaporation due to 
convection is treated by multiplying the quiescent gas evaporation rate by a 
function of drop Reynolds number (Bellan and Harstad, 1987c). The radial 
variation of the drop temperature Is estimated by solving the heat conduction 
equation within the drop and coupling It through boundary conditions to the gas 
phase (Bellan and Harstad, 1987c). Thus, for one drop 
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dm<l 
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-4irR 2 a e 
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p sat ( T s) " p Fs 


^conv 


(i) 


where f conv (Re) is the convective correction and Pp s is the fuel vapor partial 
pressure. The energy equation combined with the boundary condition yields 


_1_ (4* pt C i / T, r 2 dr) - 4*R 2 Q 0 + ^ [ C» T s + L (T s ) 
dt o y dt L 


( 2 ) 


where Q g is the heat conduction flux from the gas, Ag aTg/ar, at r = R. Note 
also that 

- 4*R 2 Pi ^ (3) 

dt dt 

The gas temperature inside the sphere of influence, along with the corresponding 
heat conduction flux, is estimated by assuming that the viscosity is a function 
of the temperature in the same manner as in Bellan and Harstad, 1988. During 
the evaporation process, the composition of the gas phase changes, causing 
concomitant viscosity and heat capacity changes. Since the present model is 
directed at providing behavioral trends, not accurate quantitative predictions, 
this effect is neglected. The ideal gas law is used to relate gas temperature 
and density taking into account that the pressure is constant. 


B. Model o f the collect ion of spheres of Influence 


Unlike the above model which was based upon differential equations, this 
model Is based upon global conservation equations. These describe 
momentum and energy conservation for both phases and their Infections. 
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Of particular interest is the description of the radial notion of the drops 
in the cluster. If r is the radial position of a drop in the cluster and R(t) 
is the cluster radius, a radial similarity parameter is f - r/R, with 
0 < « < 1. For a given drop, f is a constant parameter. Radial velocities and 

the radial electric field and electrostatic force are taken proportional to «. 
(There are additional small scale electric field variations around each drop; 
these do not effect drop motion.) On a large length scale the charge density in 
a cluster is uniform. For a uniform charge, the Integral form of Gauss' law 
gives a radial field and corresponding drop force proportional to r • c R (t). 
The drop ineitial force is proportional to d^/df, or fd a R/dt a , hence both 
forces are linear in «. Global momentum eouations for the radial and axial 
velocities are formulated in a manner similar to that of Delian and Harstad, 
1987a. The cluster boundary velocity is dR/dt, the gas expansion velocity at 
the cluster edge (f-1) is u ge and the radial slip velocity is u re - dR/dt-u ge at 
the edge of the cluster. Note that on a large length scale the evaporating 
drops form a uniform mass source density and that this is consistent with the 
similar (u re radial slip velocity. In general, the momentum equations are 
consistent with similarity with the exception of nonlinear drag and convective 
derivative terms. Average equations are formed by Integrating over parameter < 
to be consistent with the other global equations where the radial dependence 
does not appear. (This eliminates some coupling between axial and expansion 
velocities due to convective derivative terns.) This means that a tendency of 
nonlinear drag to destroy self-similar radial motion is being neglected. The 
radial drag Is proportional to [Cg (Ura a * F a J re a ) l ^ 3 l ^ u re* FUminatl g 
5 factor in the forces, the expression In the square brackets Is being averaged. 
In the small Reynolds number Stokes flow 11.1t, this bracket Is Independent of 
{i for larger Reynolds niabers the variation Is equivalent to a power of f 
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less than unity (Bel lan and Harstad, 1987b; True, 1983). This averaged drag is 
considered adequate for the purposes of the current model, which emphasizes 
comparative behavior. Due to use of self-similar radial velocities, the mass 

loss from a cluster is calculated in a consistent manner, without the use of an 

nwhiun 1988. where 

assumed "trapping factor'- as was uo».e ... ud-n na.-.a-, 

interpolation was performed from a strictly steady gas phase (maximum loss) 

limit to the null loss limit. This also means that the mass loss no longer 

depends on dR/dt, which was only approximated in Bellan and Harstad, 1988. (The 

models produce R(t) directly, hence its derivative is not directly calculated.) 

In a manner similar to Bellan and Harstad, 1988, two forms of the Nusselt number 

are used in the model. In the first, turbulence at the edge of a dense cluster 

builds up with time as the cluster porosity (axial slip velocity) collapses, 

whereas in the second, the turbulence is initially present. 


The global conservation equations are as follows: 


Conservation, 
ft?* Ifygl ± air). 

dm g . -N + 4*R 3 pn a u re + -rR 2 (pg. - Pga) u ra 

dt dt 

Ml 


- 4«R 2 p ga (1 - Ypa) “re + 
dt 

4«R ^ f ^ 1 Nuc (Y Fa - Y F J ♦ 
Pr l iki J 
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rR 2 


where 


(1 - Vp.) - Pga (1 - Y Fa ) j ”ra 

(5) 

1 dR / dt - Ugg 

(6) 

U d - “g 

(7) 


are slip velocities. 

Momentum Conservation 

Drops (cluster) 


m dun 1 _ 

I'd _ - - _ 7g 

dt 2 y 


rR*C 0 u r u ra + *R 2 C 0c ?L. u 1J!2 

(Naj + mg) 


( 8 ) 


where u r 2 ■ u ra 2 + | u re 2 , u r is the mean slip velocity magnitude, and 

pg « mg / V c , V c - 4* R 3 /3. The two terms in the right hand side of eq. 8 
represent respectively drag on an individual drop and drag on the entire cluster 
(Bel lan and Harstad, 1987a) 


Axial slip 


md 


du 


ra 


dt 


Radial slip 


md tan 

dt 


f N 5 1 dmd f Nm d + mg 1 * R2 c 

l mg J dt [ 2 V C J 


Nqd 2 r - 2 + 


J ra 


(9) 
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( 10 ) 


9 


Energy Conservation 


Gas 


where 


dE, 


- -4*R‘ NQg - N ^ Cn F T s + 
dt y dt P 

4 ’ rRw Pga c p *ga u re + *R 2 c p (pg« T gco - p ga T ga ) u ra + 

4 '" % ( £ ] Nuc C P ‘V - Tga) + 

\ *R 2 Pga u re “ge 2 ■ I N ^ u ra 2 + 


dt 


- 7g 
2 a 


wR 2 N 


C D u r u ra 2 + *R 2 C 0c Ud Ug [ Klfl d u d + "g u g ) 

l Nrad + mg J 


Eg - Hg + l (Nq d ) 2 R " 1 + 


5 4*£/ 


(11) 
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Nn>d (dR/dt) 2 + mg Ug e s 


( 12 ) 


The last two lines in the right hand side of eq. 11 are small dissipation terms 
proportional to the square of the Mach number. Likewise, the last term in E g is 
small. The derivative dR/dt is approximated by a third order backward 
difference; it appears only in terras that are nearly negligible. 

The total cluster gas mass, mg, total air mass, ra a , and total thermal 
enthalpy, H g , are given by 

Og ■ Pga V c fg 

“a M pga v c f a 

Hg - pga C p T ga V c f h 
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pi 


where f x (R, R, T s , T ga , Yp s » Yp a ) are known functions of order unity (Bel 1 an 
and Harstad, 1988). Since the pressure Is constant, the main dependency is H g - 
V c . Consequently, E g is mainly a function of R and eq. 11 may be considered as 
an equation for R. Equations 4 and 5 determine p ga and Yp a , hence also T ga . 
Equations 1 and 2 o^termine T s and Yp s . Equations 8, 9, 10 determine 
velocities. Note that the electrostatic effects appear as the first term on the 
right hand side of the radial slip momentum equation 10 for the electrostatic 
force, and the second term in equation 12 for the electrostatic field energy. 

III. Results 

All results presented here pertain to n-decane drops evaporating in an 
unvitiated ambient air. The ambient air temperature is 1000 K, the initial 
drop temperature is 350 K and the pressure is one atmosphere. The Initial drop 
velocity is 5 m/s and the initial drop radius is 20 pm. The value of other 
parameters used in the calculation are given in Bellan and Harstad, 1987c. 
Similarly to Bellan and Harstad, 1918, the form of the cluster Nusselt number 
wherein turbulence builds up with time for dense clusters will be called 
turbulence model 1; turbulence is prisent initially for turbulence model 2. 

In Figure 2, a characteristic evaporation time Is shown versus initial 
air/fuel mass ratio for a 10 cm radiis cluster with different drop charges. The 
time plotted is that when the drop radius decreases to 0.3 of its original size 
(drop mass 97.3% evaporated). The drop charge shown Is relative to the mean 
value given by Kelly, 1984: 

(<ld)max - (7-34 x 10’“ coul/cmJRo . 
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This charge value depends only on the initial drop radius. The value apparently 
applies to all hydrocarbons and appears to be the maximum charge obtainable for 
practical sprays (Kelly, 1984). Rote that for drops charged to the Kelly value, 
the Rayleigh charged drop shape instability occurs when the drop radius becomes 
somewhat smaller than 0.3 original size. The calculation Is stopped as this 
instability is encountered, avoiding modeling of the ensuing drop breakup. 
Fluid dynamic forces act in concert with drop charge In the disruption of drops 
(Cerkanowicz, 1981). For the conditions of the present calculations, this 
effect is small, typically causing a reduction in the Rayleigh limit charge by 
less than 10%. 

A spray whose initial air/fuel mass ratio is greater than or equal to 1.57 
will be termed "nondense," with a corresponding maximum drop density of 
approximately 10* drops/cm®. (This air/fuel mass ratio is 0.1 of 
stoichiometric.) The concept of denseness or cluster porosity directly 
corresponds to the relative amount of velocity slip of the ambient gas with 
respect to drops in the cluster. Comparison of prior and current calculations 
shows that for nondense clusters, the nonnegllble axial slip velocity depends 
very weakly oi, the dr*p charge, initial cluster size, and turbi'ence model. In 
contrast, clusters having slip velocities of order 10 cm/s or less are 
considered dense. Thus a cluster denseness or porosity criterion based on 
initial air/fuel mass ratio only, Independent of the above parameters, is 
expected to be valid. As Is seen In Figure 2, only dense sprays are cveatly 
affected by drop charging. The charge acts to expand the cluster into the hot 
ambient air, promoting evaporation. 

The volume expansion ratio at the characteristic evaporation time is shown 
in Figure 3. Note that for dilute sprays the electrostatic forces are too weak 
to produce any significant expansion. As was shown in Bell an and Harstad, 1988, 
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| evaporative gas cooling leads to cluster contraction for uncharged drops in 

dense clusters. The expansion effect of the radial evaporation flux is very 
J weak, being proportional to the square of the small flow-Mach-number as can be 

stsn from the terms proportional to u ge 2 in eq. 12. 

Prior and present calculations show that the evaporation time does not 
depend on the turbulence model for nondense sprays, independently of whether the 
drops are charged or uncharged. For dense sprays with significant charge (one 
quarter or more of the Kelly value), the evaporation time is slightly less with 
turbulence model 2, with a corresponding slightly greater volume expansion 
ratio. In contras*, the turbulence model can greatly affect uncharged dense 
sprays, as seen in Bellan and Harstad, 1988. This is due to the limited 

available thermal energy for evaporation in the relatively small drop 
interstitial region of dense sprays, along with comparitive isolation from the 
hot ambient. Only turbulence breaks this isolation for uncharged dense 
clusters, whereas charge- induced expansion into the ambient dominates for highly 
charged drops. For a dense spray with initial air/fuel mass ratio of 0.471, the 
behavior of the drop evaporation time relative to drop charge for the two 
turbulence models is snown in Figure 4. The time is independent of turbulence 
model for drop charges greater than half the Kelly value for this case. For 
turbulence model 1 at low drop charge, saturation before complete evaporation 

occurs. Table I presents a summary of the control parameters in various 
regimes. 

The relative amount of fuel ejected from the drop cluster as a function of 
initial cluster size is presented in Figure 5 for an initial air/fuel mass ratio 
of 1.57, the dividing value for dense/nondense sprays under the initial 
conditions given in the first paragraph of this section. Curves are given *or 
null and maximum drop charge for the two turbulence models. The fuel loss ratio 
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is an important quantity since for dense clusters ignition outside of the drop 
cluster is expected (Bellan and Harstad, 1987b), with only the ejected fuel 
participating in ignition. The fuel loss ratio depends strongly on cluster 
size, fairly strongly on the turbulence model, and least, but still 
significantly, on drop charge. Unlike the fuel loss, the calculations show that 
the total mass loss is nearly independent of the turbulence model, with a very 
weak dependence for small clusters only. Also, for charged drop clusters, the 
relative amount of mass entrainment of the ambient upon expansion increases with 
decreasing cluster size. The larger fuel loss and greater entrainment for 
smaller clusters may be attributed to their larger surface to volume ratio. 
(The curve in Figure 5 for turbulence model 1 with no charge follows an 
approximate inverse radius dependency). 

Some limited results of charge effects on ignitiun are given next. For 
nearly dense drop clusters, with initial air/fuel mass ratio greater than or 
nearly equal to 1.57 but less than 5, the penetration of the ambient is small 
enough to allow for the possibility of a diffusive ignition criterion (Bellan 
and Harstad, 1987a, 1987b). In this range the iv,.iition time decreases with 
increasing drop charge (drop radius at ignition increases). These effects are 
at most a factor of two, and largest at the higher cluster densities. Moreover, 
the ignition time increases as drop number density increases; dense clusters in 
general do not ignite before nearly complete drop evaporation (see also Bellan 
and Harstad, 1937b). Since the ignition time decrease* with increasing drop 
charge the salutory effect of drop charge in reducing evaporation time may be 
lessened by the ignition and burning of larger drops in nearly dense clusters, 
which could produce an increase in unwanted byproducts such as cenospheres for 
multicomponent fuels which are not Diesel fuels. 
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Mow the above results can ba used to infer the effect of electrostatic 
charging upon soot control by using the experimental results of Sangiovanni and 


Liscinsky, 1984. 

Despite the fact that the present calculations are made for n-decane and 
the results of Sangiovanni and Liscinsky, 1984, are obtained for n-hexadecane, 
qualitative comparisons can still be made because Figure 2 of Sangiovanni and 
Liscinsky, 1284, shews that although at constant initial droplet spacing the 

soot emission index decreases with decreasing aromatic content of the fuel, the 
general trend of the variation of the soot emission index with drop spacing is 
the same, independent of the fuel. The quantitative reduction in the value of 
the soot emission index with increasing drop spacing found by Sangiovanni and 
Liscinsky, 1984 should be considered as a minimum reduction since theirs is a 
single-droplet-stream experiment and thus the impact of increasing the spacing 
of lateral adjacent drops cannot be assessed. With this restriction in mind, 
the results of Fig. 2 of Sangiovanni and Liscinsky, 1984, are directly 

comparable to our results because the relative drop spacing here is in the same 
range of values as in that work. Here R 2 0 measures the same distance as A/d 0 in 
Sangiovanni arid Liscinsky, 1984, which was varied there from 5 to 20. In the 
present calcuations for - 0.471, T°g a » 1000 K and Y°Fva ■ 0 the cluster is 
characterized by n° - 3.25 x 10 4 cm* 8 and R 2 0 * 8.8; and for <t>° - 4.71 at the 
same conditions, we obtain n° - 2.98 x 10 8 cnT 8 and R 2 0 - 19.6. The first case 
corresponds to a dense cluster of drops where interactions are strong and the 
second case corresponds to a dilute cluster (although a rich mixture) where 
interactions are weak (see Fig. 2 for the uncharged limit). Although the 

difference between the ambient temperature used in this study and that of 
Sangiovanni and Liscinsky, 1984, might somewhat influence the results 

quantitately but not qualitatively, the conclusions should be independent of the 



mass percent of ambient oxygen (in the low range) as shown in Fig. 3 of 
Sangiovanni and Liscinsky, 1984. Since all mass fractions must add to unity, we 
conclude that the results should also be independent of V°Fva in the 
corresponding range. Previous results obtained for uncharged drops (Bel lan and 
Harstad, 1988) show that the radius of the sphere of influence can either 
decrease for very dense clusters of evaporating drops (due to cooling effects), 
or remain nearly constant in the case of dilute clusters. In contrast, when 
charging, Fig. 3 shows that the radius of the sphere of influence, and thus that 
of the cluster, increases; for example when * 0.471 and the charge ratio is 
unity (the largest charge feasible), R °2 » 8.8 and R 2 * 17.2 at Rj * 0.3 (the 
radius near Rayleigh instability). There is a continuous cluster expansion 
(monotonic variation) starting from the initial condition and therefore R 2 is 
larger than R 2 0 ; thus, the results of Sangiovanni and Liscinsky, 1984, can be 
used here uncontroyersially. According to Fig. 2 of Sangiovanni and Liscinsky, 
1984, an increase in R 2 from 9 to 17 (compare with the above 8.8 to 17.2) could 
result in a substantial decrease in the soot emission index. Since Fig. 3 shows 
that at constant charge the volume ratio increases with decreasing 4 >° (and thus 
R 2 0 ), electrostatic spray dispersion yields more dramatic results at constant 
charge as the cluster is denser. 
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| IV. Summary and Conclusions 

The above results shuw that, unlike the case of uncharged drops in dense 
] clusters (Bellan and Harstad, 1988), when the drops are charged, turbulence is 

not a control parameter for dense clusters of drops; instead electrostatic 
effects dominate. For dilute clusters of drops it is found that similarly to 
the uncharged drops (Bellan and Harstad, 1988), turbulence effects are 
unimportant. Significantly, the results above show that the evaporation time of 
dilute clusters of drops cannot be changed by charging the drops in contrast to 
the dense clusters of drops where the evaporation time can be considerably 
shortened by this procedure. Moreover, it is also found that the fuel vapor 
loss from the cluster increases as turbulence increases and as cluster size 
and/or charge decreases. 

It was also shown by comparing with existing experimental results that 
when charging the drops with the maximum empirical limit obtained by Kelly, 
1984, substantial decrease in the soot emission index can be achieved for dense 
clusters of drops. These results are encouraging and further studies of 
electrostatic spray dispersion should be made in order to ascertain its full 
range of posibil ities. 
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Symbols 


Nomenclature 



Y 




a 

^D» ^Oc 
c i 

c p> Cpp 
1 1 

ng 

k G 

L 


m 


n 

N 


Nuc 

P 

Pr 

<ld 

r 

R 

R 

R 1 

r 2 

Re 

t 

T 


J ge 

J ra 


J re 


Wf 


*Fa 

Yf« 


radius of the sphere of influence of each drop 

drag coefficients for drop, cluster 

liquid specific heat 

specific heat of gas, fuel 

total enthalpy of cluster 

universal gas constant 

latent heat of evaporation of fuel 

mass 

drop number density 

number of drops 

Nusselt number of cluster 

pressure 

Prandtl number 

drop charge 

radial coordinate 

drop radius 

cluster radius 

R/r<° 

a°/R° 

Reynolds number of a drop 
time 

temperature 

velocity 

gas expansion velocity at cluster edge 

axial slip velocity 

expansion slip velocity 

cluster volume 

fuel molar weight 

interstitial fuel mass fraction 

ambient fuel mass fraction 
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-A 




J 


-2-l-A r JL- 


/ ”, 


J , 


J 


J 


l- 


aC % 


j 

o 


Greek Symbols 


o 

4 




evaporation accomodation coefficient 

permittivity of free space 

gas thermal conductivity 

ambient viscosity 

iiidss density 

air fuel, mass ratio 


Subscripts 

a 

d 

9 

ga 

g« 

i 

r 

s 


drop interstitial region 
drop 

gas in cluster 
interstitial region gas 
ambient gas 
liquid 

relative or slip velocity 
drop surface 


Superscripts 

o 


initial condition 
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TABLE 1 


Control Parameters for the Evaporation Time 


Control Parameter 
for the Evaporation 
Time 

*• - ~ i- i j ( > \ 

uncntryeu urup**-* 

1 *JI 

Liiai ycu 

drops 

dense 

non -dense 

dense 

non dense 

Turbulence 

important 

not important 

not 

important 

not 

important 

Electrostatic 
charge when 

— 

— 

important 

not 

important 

dd > ifld-j-max 

4 






(a) Bel lan and Harstad, 1988. 
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Figure Captions 


Fig. 1 Moving Spherical Cluster. 

Fig. 2. Characteristic evaporation time vs initial air/fuel mass ratio. Drop 

charge is relative to maxima®, see text. 

Fig. 3. Volume expansion ratio vs Initial air/fuel mass ratio. Drop charaa is 
relative to maximum. 

Fig. 4. Characteristic evaporation time vs ratio of drop charge relative to 
maximum. Spray cluster Is dense, see Bel lan and Harstad, 1988, for 
turbulence models. 

Fig. 5 Relative fuel ejection from cluster vs initial cluster radius. 

Cluster density is at the boundary of the near-dense/nondense regimes. 
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